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Abstract 

For the first time, the tilted axis cranking covariant density functional theory with pairing 
correlations has been formulated and implemented in a fully self-consistent and microscopic way 
to investigate the evolution of the spin axis and the pairing effects in rotating triaxial nuclei. The 
measured energy spectrum and transition probabilities for the ^^®Nd yrast band are reproduced well 
without any ad hoc renormalization factors when pairing effects are taken into account. A transition 
from collective to chiral rotation has been demonstrated. It is found that pairing correlations 
introduce additional admixtures in the single-particle orbitals, and, thus, influence the structure of 
tilted axis rotating nuclei by reducing the magnitude of the proton and neutron angular momenta 
while merging their direction. 
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I. INTRODUCTION 


Similar to rotational bands observed in molecules, the most common collective excita¬ 
tion in nuclei corresponds to a rotation about the principal axis of the density distribution 
with the larger moment of inertia. The rotation is collective since a large fraction of the 
angular momentum is generated through small contributions from many nucleons. It is 
well known that nuclear deformation and superfluidity play vital roles in generating angu¬ 
lar momentum [l|. The substantial deformation of the overall density distribution specihes 
the orientation of a nucleus and, thus, the rotational degree of freedom. Meanwhile, its 
superfluid behavior is required by the fact that the observed collective moment of inertia is 
usually much smaller than the rigid-body estimate. 


Unlike molecules, nuclei can rotate about an axis tilted with respect to the principal axes 

n 

of the density distribution |2[ due to the fact that a nucleus is composed of nucleons carrying 


a quantized amount of angular momentum. This is the so-called tilted axis ro 


ation which 
4|. The tilt 


was hrst proposed within the mean-held tilted axis cranking (TAG) approach 
of the rotational axis is closely related to the interplay between collective and single-nucleon 
motions. Therefore, the two elements, deformation and superhuidity, are not only crucial in 
generating the magnitude of the spin but also its orientation. 

A variety of discrete symmetries for rotating nuclei can be obtained by combining the 
deformation and the spin orientation, and this gives rise to a variety of new phenomena. 
Among the latter figure, the magnet.c aud autimaguetie rotation in nearly spher.eal nude. Q 
a, 1^, the high-A' bands giving rise to K isomerism [7[ in axially deformed nuclei, etc. For 
triaxial nuclei, specihcally, the rotational axis may lie outside the three principal planes of the 
ellipsoidal shape. This forms the so-called aplanar rotation and causes nuclear chirality js], a 
mode that has attracted signihcant attention due to its importance on the subatomic physics 
scale { 2 , 9|. For axially deformed nuclei, the rotational axis always lies in the plane dehned 
by the symmetry axis and the one perpendicular to it. In contrast, planar triaxial solutions 
can be found in the three distinct principal planes. This also leads to many new interesting 
modes such as chiral vibration ly, lll| and transverse wobbling 12l-ll4|. Furthermore, a 
rotating triaxial nucleus allows more degrees of freedom for the evolution of the spin axis. 
This makes it all the more interesting to investigate, for example, a transition from planar 
to aplanar rotation, i.e., from chiral vibration to static chirality jlOl . Illl . Il5|. and a transition 
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from one planar rotation to another; this has not been studied in any detail so far. 
In contrast to the well-known impact of pairing on principal axis rotation 


16| . its influ¬ 


ence on a tilted axis rotor is still far from being understood. In particular, because it affects 
both the collective and valence nucleon motions, it should be expected that pairing corre¬ 
lations may reorient the spin axis in tilted rotation. So far, most calculations are based on 
single-particle potentials combined with the pairing plus quadrupole-quadrupole model jd]. 
Therefore, self-consistent methods based on more realistic two-body interactions are required 
for a more fundamental investigation, including all important effects such as core polariza¬ 


tion and nuclear currents 


17 


19j . Such calculations are more cha 


the framework of both relativistic 


20 


22| and nonrelativistic 


17 


lenging, but are feasible in 


23 


24j | density functional 


theories (DFTs). However, pairing correlations have not been taken into account in any 
of these studies. For many years, the computed bandhead energies and magnetic dipole 
transition probabilities B{M1) had to be renormalized by ad hoc factors to reproduce the 
data 3, 22, 251, leading to a long-standing question about whether the inclusion of pairing 
correlations would improve the agreement between data and calculations. 

The focus of the present research is two-fold: 1) the transition of the spin axis from 
one principal plane to another in a triaxial nucleus; 2) the impact of pairing correlations 
in a tilted axis rotor. This article presents the hrst tilted axis cranking covariant DFT 
with pairing correlations. The spin axis evolution and the impact of pairing for a rotating 
triaxial nucleus have been investigated in a fully self-consistent microscopic way for the hrst 


time. Covariant DFT 


26|-|28| consistently treats the spin degrees of freedom, includes the 


complex interplay between the large Lorentz scalar and vector self-energies induced at the 
QCD level 29(]. Moreover, the nuclear currents which are essential for rotating nuclei, are 
provided naturally from the spatial parts of the vector self-energies. 


II. THEORETICAL FRAMEWORK 

For a unihed and self-consistent treatment of the mean helds and pairing correlations, one 

n 

has to solve the fully relativistic Hartree-Bogoliubov (RHB) problem [30| in the framework 
of superhuid covariant DFT. The RHB model contains two average potentials: the mean 
helds S{r) and V^{r) which include all the long range particle-hole {ph) correlations, and 
a pairing held A(r) which sums up the particle-particle {pp) correlations. In the TAG 
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model, these potentials are deformed and the calculations are carried out in the intrinsic 
frame rotating with a constant angular velocity vector uj pointing in a direction which is 
not parallel to one of the principal axes of the density distribution: 


( h — u: ■ J A 

I -A* -h* + UJ ■ J* 



Here h = — X is the single-nucleon Dirac Hamiltonian 


( 1 ) 


ho = OL ■ (p - V) + /3{m S') 


( 2 ) 


minus the chemical potential A, and J is the total angular momentum of the nucleon spinors. 
Uk and 14 are the quasiparticle Dirac spinors and Ek denotes the quasiparticle (qp) energies. 
The mean helds S and as well as the pairing held A are connected in a self-consistent 
way to the densities and currents as well as to the pairing tensor distributions. The iterative 
solution of these equations yields expectation values of the angular momentum, total ener¬ 
gies, quadrupole moments, transition probabilities [Ml, etc. The magnitude of the angular 
velocity lo is connected to the angular momentum quantum number I by the semiclassical 
relation (J) • (J) = /(/ -|- 1). 


The observed yrast band in the odd-H nucleus ^^^Nd 3^ is investigated in the present 
work. The ground band is associated with the i'hii /2 one quasi-neutron conhguration [33 1. 
However, above / = 29/2h, this conhguration is further coupled to two aligned hii /2 protons, 
thereby resulting in the 3-qp conhguration i/hn / n. This 3-qp band and its partner 

have been interpreted as a pair of chiral bands [ul, IlSf, a property supported by lifetime 
10|. The present self-consistent investigation includes both the 1-qp and 


measurements 


3-qp bands, and the evolution of the rotational axis will be analyzed. The point-coupling 


energy density functional PC-PKl 3^ is adopted in the ph channel, and a monopole pairing 
force with constant strength, for neutrons = 0.12 MeV fm^ and for protons Gp = 0.13 
MeV fm^, determined from the odd-even mass diherences is used in the pp channel. The 
calculations are free of additional parameters. 

In this work, we allow only rotations around an axis in the {x, z) plane. Equation ([T]) 
is solved in a three-dimensional Cartesian harmonic oscillator basis 35|] with N =10 major 
shells. Parity is the only good quantum number, and the space of the Hamiltonian matrix 


is, thus, twice as large as for the principal axis cranking RHB theory 


36|. Therefore, 
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parallel computations have been implemented to reduce the required computational time. 
Another difficulty is that one has to trace and block the right qp orbitals to keep the 
multi-qp conhguration unchanged while solving Eq. ([T]) iteratively with different A and oj 
values. To achieve this, we hrst dehne for each single qp state the normalized wave functions 
i’k = Uk/^{UkVUk and rk = Vk/^/iWVk. Then, we search at each iteration for the 
largest overlap 

where k[ are determined by running over all the single qp states obtained in the previous 
step while maximizing the overlap. Finally, the state k would be blocked only if Oi is larger 
(smaller) than O 2 and the state k[ {k' 2 ) was blocked (unblocked). 


III. RESULTS AND DISCUSSION 

The present self-consistent calculation for ^^^Nd leads automatically for the bandhead to a 
state with considerable triaxial deformation (7 ~ 22°) associated with the 1-qp conhguration 
vhiij 2 - Along the band, the tilt angle 9 of the rotational axis is determined in a self-consistent 
way by minimizing the total Routhians (Fig. [1]). Here, the tilt angle is dehned as the angle 
between the rotational and the long axis, and the positive (negative) value denotes a tilt 
towards the short (intermediate) axis. Figured] clearly indicates that the 1-qp state always 
favors rotating along an axis in the long-intermediate {l-i) plane, and the rotation axis tilts 
more and more appreciably toward the i axis with increasing frequency. Moreover, the 3-qp 
state, rotating along an axis in the long-short (Z-s) plane, becomes lower in energy than the 
1-qp conhguration at high frequency. This is consistent with the fact that the observed 3-qp 
band becomes yrast at high angular momentum. The pairing ehects for the 3-qp states are 
suppressed signihcantly by the two aligned hii /2 protons. In contrast, the 1-qp states are 
lifted notably in energy by excluding the pairing correlations, and, thus, they are always 
located much higher than the corresponding 3-qp states. Therefore, it is concluded that 
one can determine the band crossing frequency more accurately by including the pairing 
correlations. 

In addition to determining the tilt angle, the calculated energy spectrum and the angular 
momenta could be calculated as well: these are compared with the data in Fig. [2| with 
the upper panel indicating that the experimental rotational excitation energies for both the 
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FIG. 1. (Color online) Total Routhians for the 1-qp and 3-qp configurations as functions of the tilt 
angle 9 at the rotational frequencies fvjo = 0.2 (top), 0.3 (center), 0.4 (bottom) MeV. The circles 
and triangles represent the local minimum points for the 1-qp and 3-qp configurations, respectively. 

1-qp (lower spin part) and 3-qp bands (higher spin part) are reprodnced well by the present 
self-consistent calculations with pairing. In particular, the inclusion of pairing leads to the 
correct energy difference between the 1-qp and 3-qp conhgurations. Therefore, the need for 
an artihcial renormalization of the bandhead energies mentioned above has been eliminated 
by including pairing. 

In the lower panel of Fig. |2l it is seen that the calculated total angular momenta also agree 
well with the data. A backbending happens in the region I = 10-14h, where the angular 
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FIG. 2. (Color online) Rotational energy (upper panel) and rotational frequency (lower panel) as 


functions of the angular momentum in comparison with the data of Ref. 


15l | (solid dots). Here, 


the excitation energies are the relative energy differences with respect to the ground state. 


momentum increases while the rotational frequency drops drastical 
such a phenomenon is beyond the scope of a cranking calculation 


It is wellknown that 


37( 1 ■ and, consequently, 


the calculated results are omitted. It should be expected that the moment of inertia J/o; is 
reduced by pairing in the present tilted axis calculations. However, it is surprising that the 
reduction for the 3-qp band is only visible at high frequency, where the pairing effects tend 
to be hindered by the Coriolis term. Indeed, here the proton pairing gaps are very small 
due to the two quasi-protons, while the neutron ones decreases from 1.25 MeV {hjjj = 0.2 
MeV) to 1.14 MeV {hw = 0.6 MeV). 


To understand this distinctive feature, the neutron and proton angular momentum vectors 
are provided in Fig. |3]for both the 1-qp and 3-qp bands. For the former band, the neutron 
angular momentum aligns along the I axis and the proton one essentially vanishes at fkn = 0.1 
MeV. Along the band, the neutron angular momentum keeps its projection on the I axis 
nearly constant, which reflects the contribution of the one unpaired neutron hole in the hii /2 
shell, but at the same time the proton angular momentum increases remarkably along the i 
axis due to a coherent collective rotation of many nucleons. 

For the 3-qp band, the neutron angular momentum aligns mainly along the I axis due 
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FIG. 3. (Color online) Neutron J^, and proton angular momentum vectors for both the 1-qp 
(left) and 3-qp (right) bands. Different types of arrows denote the results at different rotational 
frequencies. Note that for the 1-qp band is negligible at hw = 0.1 MeV. 

to the hii /2 hole, while the proton one aligns mainly along the s axis, due to the two hii /2 
particles involved. As the frequency increases, the neutron and proton angular momenta 
align toward each other and generate larger total angular momentum with the direction 
nearly unchanged. This situation is reminiscent of the “shears mechanism” in a magnetic 
rotation band. However, due to the considerable triaxiality (7 ~ 22°) involved, it has been 
shown in Ref. that a chiral vibration, resulting from rapid conversion between the left- 
handed and right-handed configurations, has been realized in the 3-qp band. Therefore, a 
transition from collective to chiral rotation is observed here. 

In contrast to principal axis rotation, pairing effects on the angular momentum here 
result in two competing effects. On the one hand, the magnitudes of both proton and 
neutron angular momenta are reduced by pairing. On the other, however, pairing tends to 
reduce the angle between the proton and neutron angular momenta (by up to 20 %) and 
in this way increases the total spin. Therefore, the impact of pairing on the total angular 
momentum can be, in some cases, counteracted as is the case for the lower spin part of the 
3-qp band in Fig.O All in all, pairing introduces superfluidity in rotational states and, here, 
it expedites the closing of the proton and neutron angular momenta. 

To trace the microscopic reason for the pairing effects, it would be quite helpful to use the 


Bogoliubov transformation to transform the qp basis to the canonical (particle) basis 16| . 
as in a microscopic picture, the angular momentum comes from all the individual particles. 
Since the pairing effects on the protons are significantly blocked in the 3-qp band, we show 
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FIG. 4. (Color online) Angular momentum alignments for the 1-qp and 3-qp configurations on the 
z axis (upper panels) and x axis (lower panels) of the neutron holes in the hii /2 shell (left) and 
the ( 57 / 2 c? 5 / 2 <i 3 / 2 'S 1 / 2 ) shell (right). 


in Fig. IHonly the angular momentum contributions of the neutron holes while noting that 
the proton ones lead to a similar conclusion. 

For the 75 neutrons in there are contributions from only the seven neutron holes 

with respect to the closed N = 82 shell; these include three negative-parity ones in the 
hii /2 shell and four positive-parity ones in the {g 7 / 2 d 5 / 2 dz/ 2 Si/ 2 ) shell with low j values. One 
unpaired hole always occupies the hii /2 shell and, thus, the alignment along the 2 ; axis, i.e., 
I axis, is almost constant. As uj becomes larger, the increase in angular momentum is gener¬ 
ated mostly along the x axis, i.e., i axis for the 1-qp band and s axis for the 3-qp one, through 


mixing of orbitals with large jx components. This is similar to the mec 


lanism pointed out 
3l| . Pairing always 


in the previous self-consistent investigations without pairing of Refs. 3, 
reduces the alignment of the low-j orbitals in the gds shell, and this is connected with the 
fact that it reduces the collective moment of inertia by allowing partial occupation of the 
single-particle orbitals. Moreover, it is interesting to note that the alignment Jx for the 
three i^(hii/ 2 ) holes rises with the inclusion of pairing, which indicates that additional ad¬ 
mixtures are also introduced to the valence particles (holes) orbitals. Therefore, it is clear 
that the self-consistent nucleonic Cooper-pair dynamics of pairing correlations influences the 
single-particle orbitals and their occupation probabilities, i.e., allows additional mixing in 
the single-particle orbitals, and, thus, influences the generation of the nuclear spin. More- 
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over, due to the strong Coriolis term at high frequency, pairing effects become weaker and, 
thus, the increment of the alignment Jx becomes smaller as well. 
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FIG. 5. (Color online) Calculated B{M1) (upper panel) and B{E2) (lower panel) values as func¬ 
tions of the angular momentum in comparison with data [l^ (solid dots). 


Finally, the calculated transition probabilities are given in Fig. [5] in comparison with 
the available data lOj]. Good agreement is achieved after pairing correlations are included, 
especially for the 3-qp band. The B{M1) values are derived from the relativistic expression 


of the electromagnetic current operator, which includes 


rents 


2l[ |. Since the hrst relativistic TAG calculation 


20 


Doth the Dirac and anomalous cur- 
, an artihcial factor of 0.3 has been 


used for many years to attenuate the oversized values of Here, however, it is found 

that this factor is not needed as an increase in B{M1) probabilities is counteracted mostly 
by pairing, which reduces the transverse magnetic moment by merging the directions of the 
proton and neutron angular momenta, as shown in Fig. [3l In addition, it should be noted 
that the nuclear deformation is almost unchanged by including pairing. The slight rise of 
the B{E2) values is mainly due to the fact that pairing moves the rotational axis slightly 
away from the I axis, i.e., there is more susceptibility to rotational alignment. 


IV. SUMMARY 

In summary, the hrst tilted axis cranking covariant DFT with pairing correlations has 
been formulated and the evolution of the spin axis in the yrast band of ^^^Nd has been 
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investigated in a fully self-consistent microscopic way. The longstanding question of how 
pairing correlations influence the structure of tilted axis rotating nuclei has been addressed. 
The present work shows that the experimental energy spectrum and the transition proba¬ 
bilities are well reproduced when pairing effects are taken into account. In particular, the 
artihcial renormalization for both the bandhead energy and the B{M1) values is eliminated 
by including pairing. Moreover, it is found that the superfluidity induced by pairing allows 
additional mixing in single-particle orbitals, and influences the generation of the total spin, 
i.e., reducing the magnitude of the proton and neutron angular momenta, but expediting 
the merging of their directions. 
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